Dynamic sensitivity of photon-dressed atomic ensemble with quantum criticahty 
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We study the dynamic sensitivity of an atomic ensemble dressed by a single-mode cavity field 
(called a photon-dressed atomic ensemble), which is described by the Dicke model near the quantum 
critical point. It is shown that when an extra atom in a pure initial state passes through the cavity, 
the photon-dressed atomic ensemble will experience a quantum phase transition, showing an explicit 
sudden change in its dynamics characterized by the Loschmidt echo of this quantum critical system. 
With such dynamic sensitivity, the Dicke model can resemble to the cloud chamber for detecting a 
flying particle by the enhanced trajectory due to the classical phase transition. 
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I. INTRODUCTION 

The quantum phase transition (QPT) fH occurs at zero 
temperature when the external parameters of some in- 
teracting many-body systems change to reach the crit- 
ical values. Generally, it is associated with the ground 
state with energy level crossing and symmetry breaking 
at the critical points. Recently, it was discovered that, 
near the quantum critical point the QPT system pos- 
sesses the ultra-sensitivity in its dynamical evolution Q . 
This theoretical prediction has been demonstrated by an 
NMR experiment [3]. Similar sensitivity exists in some 
quantum systems jJ49j possessing QPT. 

In this paper, we study the dynamic sensitivity of an 
atomic ensemble in a cavity with a single-mode electro- 
magnetic field (called a photon-dressed atomic ensem- 
ble), which is described by the Dicke model [l^i- We 
assume the atoms in the Dicke model are resonant with 
the cavity field. When an extra two-level atom in large 
detuning goes through the cavity field, the frequency of 
cavity field will be shifted effectively according to the 
Stark effect so that the photon-dressed atomic ensemble 
near the QPT will be forced into its critical point. In 
this situation the dynamic evolution of the Dicke model 
becomes too sensitive in response to the passage of the 
extra atom. 

Here, this dynamic sensitivity is characterized by the 
Loschmidt echo (LE) which is intrinsically defined 
by the structure of the photon-dressed atomic ensemble. 
For a short time approximation, we prove that the LE is 
just an exponential function of the photon number vari- 
ance in the photon-dressed atomic ensemble. This finding 
means that the LE can be experimentally measured by 



detecting the photon correlation. Its sudden change may 
imply the passage of an extra atom through the cavity. 
With this reorganization we will demonstrate that such 
quantum sensitivity in the Dicke model is very similar 
to the classical sensitivity of the cloud chamber for de- 
tecting a fiying particle, which is characterized by the 
macroscopically observable trajectories enhanced by the 
classical phase transitions. 
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FIG. 1: (Color online) Schematic of a cavity field coupled 
with an atomic gas consisting of A'^ two-level atoms. An extra 
detected two-level atom S is injected into cavity field. 



This paper is organized as follows. In Sec. |TT1 we de- 
scribe the setup of the quantum critical model based on 
the Dicke model. The effective Hamiltonian is given in 
terms of the collective excitation of the atomic ensem- 
ble. Then the analytic calculation of LE (or the decoher- 
ence of the extra atom) is carried out in Sec. IIIII for the 
normal and super-radiant phases, respectively, by short- 
time approximation. In the following Sec. IIVI we plot 
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some figures to explicitly show the sensitive properties 
of the LE. In Sec. El we address the similarity between 
the dynamic sensitivity of the photon-dressed atomic en- 
semble induced by an extra atom and the classical cloud 
chamber. Finally, we draw our conclusion in Sec. I VII 
The detailed coefficients for Bogoliubov transformation 
in Sec. IV are given in the Appendix. 



formation jT3, [13 (or other elimination methods) to ob- 
tain the effective total Hamiltonian 
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II. MODEL AND HAMILTONIAN 

As showed in Fig. (TJ we consider an atomic ensemble 
confined in a gas cell coupled with a single-mode cav- 
ity field of frequency w, which is described by the an- 
nihilation (creation) operator a (a^). We use the Pauli 



matrices ai^^ = |e)jj(e| - \g)n{gV <t+' = \e)jj{g\, and 

f-'' = \9)jj{'A to describe the atomic transition of the 
jth atom with energy level spacing wq, where \e) ^ and 
\g) j are the excited and ground states of the jth atom, 
respectively. The system of the atomic ensemble cou- 
pled with the single-mode cavity field is described by the 
Dicke model (hereafter, we take h= 1), 
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Here, for small-dimension atomic gas 'S'], we have as- 
sumed that all the atoms locate near the origin point and 
interact with the cavity field with the identical coupling 
strength g^. 

An extra two-level atom S with transition operators 
CTj; , cr_|_ , and a_ couples to the original single- mode cavity 
field with Hamiltonian 



(2) 



where we have made a rotating wave approximation. 
Similarly, Ws is the transition frequency between the 
ground state \g) and excited state |e) of the atom S; 
Qs is the corresponding coupling strength. 

It has been shown that the QPT will occur in the sys- 
tem described by Dicke Hamiltonian ([1} 5] , since it keeps 
Hermitian only for a small coupling strength go. But 
it is only a model to display QPT in quantum optical 
system. Actually it could not happen for the realistic 
atomic, molecular, and optical (AMO) system if the un- 
reasonably ignored two-photon term is included ^T5]. 
To focus on our main idea in the work, we only regard 
the Dicke system as a simplified model. We would like to 
point out that many authors have recognized this prob- 
lem, but there still exist many explorations by using this 
simplified model [isj . 

If the atom S is far-off-resonant with the cavity field, 
that is, the detuning (= — w) is much larger than 
the corresponding coupling strength gg, i.e., |As| ^ g^, 
then one can use the so-called Frohlich-Nakajima trans- 



where 6 = g1/^s and g = govN- We note that the 
Frohlich-Nakajima transformation is equivalent to the 
approach based on the adiabatical elimination. 

The Hilbert space of N two-level atoms is spanned by 
2^ basis states. In the current case all the atoms have 
the same free frequencies and coupling constants with 
the cavity field, we can consider these atoms being iden- 
tical. Then the Hilbert space is reduced into a subspace 
of {2N+1) dimension. In this subspace, Hamiltonian (|3]) 
is simplified by introducing the collective atomic opera- 
tors 



J± 



N 



1 ^ 
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which obey the following angular momentum commuta- 
tion relations. 



[Jz, J±] — ^J±, [J+, J-] ~ "^Jz 



(5) 



The collective atomic operator Jz denotes the collective 
population of the atomic gas and J± represents the col- 
lective transitions. 

In terms of the above angular momentum operators, 
Hamiltonian ([3]) is written as 

Hcs = {uj + 6az)a''a+^{uJs + 5)az 

+ojoJz + ^{a^ + a){J+ + J^), (6) 
V N 

which is further reduced to 



+g{a^ + a ) ib'^^l-h%/N + h.c\ 



(7) 



(up to constant terms) through making use of the 
Holstein-Primakoff [iBl transformation, which represents 
the angular momentum operators in terms of a single 
bosonic mode as follows: 



J+ 
J- 

Jz 



bWN -b'fb, 

b'<b- 



b'lbb, 



(8) 



To see more explicitly the dynamic sensitivity of the 
photon-dressed atomic ensemble in response to the extra 
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atom, corresponding to different state of the extra atom, 
the effective Hamihonian in Eq. ([7|) reads 

Hcs = \g) {g\ ®Hg + \e) (e| ® (9) 

with 

Hg = ujga^a + LOob^b + gi^a' + a) 

X (^b^^Jl-b^b/N + h.c^ , (10) 

He = Wefl^a + a;o&^6 + .g(a^ + a) 

X (^b^^l-b^b/N + h.c}j , (11) 

where cje = uj + 6 and tUg = lu — S. Note that in the 
derivation of the above Hamihonians PH)) and (|lip . we 
have discarded some constant terms. 



that the dynamic of the photon-dressed atomic ensemble 
is sensitive to the state of the extra atom. When the 
photon-dressed atomic ensemble stays in the vicinity of 
the QPT, the effect of QPT must impose on the state 
of the extra atom with some enhancement fashion, like 
the results in Ref. Q. This motivates us to study the 
quantum decoherence of the extra atom near the critical 
point of the photon-dressed atomic ensemble, which can 
also reflect the dynamic sensitivity of the photon-dressed 
atomic ensemble. 

By tracing over the degree of freedom of the photon- 
dressed atomic ensemble in evolution state ([T3|) . the re- 
duced density matrix ps(t) = Tra,f,{|^' (t)) (* (t)|} of the 
detected atom is obtained as 

Ps{t) = lap \g) {g\ + |e) (e| + (Da* (3 |e) {g\ + h.c), 

(15) 

where we have introduced the decoherence factor 



III. QUANTUM CRITICAL EFFECT 

Before the extra atom S is sent into the cavity, the 
photon-dressed atomic ensemble (including the cavity 
field and the atomic gas) is described by the Dicke Hamil- 
tonian 

Hg — ujo'a + cjob^'b + g{a^ + a) 

X (^b''^Jl-b^/N + h.c^ . (12) 

Comparing Eqs. and pT|) with Eq. p^ . we find, as a 
result of the injection of the atom S, only the frequency of 
the optical field changes by a small shift S in the dynamic 
evolution of the photon-dressed atomic ensemble. 

The photon-dressed atomic ensemble is initially pre- 
pared in the ground state \G) of Hamiltonian (|12l) and the 
extra atom S* in a superposed state a\g)+l3 |e) , where the 
normalization condition requires |ap -I- |/3p = 1. When 
the extra atom S interacts dispersively with the cavity 
field, the total system is governed by Hamiltonians pO)) 
and pT|) corresponding to the extra atom S in states \g) 
and |e), respectively. Then at time t the state of the total 
system becomes an entanglement one. 
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«*(a|5)+/3|e))® |G) 
= a\g)®e~'"^'\G)+P\e)®e-'"''\G) 
= a\g) ®\Gg{t)) + p\e) ®\GS)) , (13) 



where we have defined 



-iH„t 



|G), |Ge(i))^e-^=*|G). (14) 



The generation of the above entanglement is due to the 
conditional dynamics of the total system. This is to say, 
corresponding to the detected atom prepared in states 
\g) and |e), the evolution of the photon-dressed atomic 
ensemble will be governed by the Hamiltonians Hg and 
-ffe, respectively. The central task of this paper is to show 



D[t) = (G| cxp{iHgt) exp {-iHet) |G) 



(16) 



Alternatively, we can investigate the decoherence of the 
extra atom by examining the so-called LE 



(17) 



defined for the dynamic sensitivity of the photon-dressed 
atomic ensemble. For a short time t, the LE can be 
approximated as 



(G| e-^'^^'"^" |G) 



(18) 



The straightforward calculation can give 

Lit) ss exp (-47(5^^2^ . (19) 
Here, we have introduced the photon number variance 



7^((ata)^)-(atar 



(20) 



and the average (•) is taken for the ground state |G). 

We point out that, up to the second order of time t, 
the decay rate of the LE depends not only on t^, but also 
on the photon number variance 7. It is well known that 
the photon-dressed atomic ensemble described by Dicke 
Hamiltonian p2p transits from the normal phase to the 
super-radiant one with the increase in the parameter g 
from that less than the critical value gc = wwq / 2 to that 
larger than gc- Going across the phase transition point, 
the ground state of the photon-dressed atomic ensemble 
experiences a complex change. We can predict that the 
photon number variance 7 of the ground state will exhibit 
some special features at the critical point. 

According to Eq. p3p . we can imagine that the quan- 
tum criticality of the photon-dressed atomic ensemble 
can display which single state \g) or |e) that the extra 
atom stays. When L{t) approaches zero, the photon- 
dressed atomic ensemble is forced into two orthogonal 
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states \Gg{t)) and \Ge{t)), and thus it behaves as a mea- 
surement apparatus to detect the state of the extra atom. 
In this case, its measurement on the atom will induce the 
decoherence of the extra atom. 

In what follows, we will calculate the photon number 
variance 7 of the photon-dressed atomic ensemble in two 
different phases, that is, the normal phase and the super- 
radiant phase. 



A. Dynamic sensitivity in normal phase 

In this subsection, we explicitly calculate 7 to investi- 
gate the properties of the LE when the photon-dressed 
atomic ensemble is within the normal phase. In the 
case of low excitations at thermodynamic limit N ^ 00, 
Hamiltonian (|12p becomes 

Hg = uja''a + iUob^b + g{a^ + a){b^ + b) (21) 



for -^/l — b%/N « 1, which is typical to describe two- 
coupled harmonic oscillators. It is well known that 
Hamiltonian (|2ip becomes non-Hermitian in the over- 
strong coupling region g > gc, namely, the Hamiltonian 
possesses imaginary eigenvalues [17^. This means effec- 
tive Hamiltonian (|2ip is ill-defined for g > g^. Therefore, 
we now restrict the Hamiltonian within the so-called nor- 
mal phase region g < gc- Correspondingly, this limited 
Hamiltonian (PT|) describes the normal phase of the Dicke 
model. 

In the normal phase, Hamiltonian (PT|) can be diago- 
nalized as 



Hg = ijJAA'^A + ijjBB'^B 



(22) 



by introducing the polariton operators A (A^) and B 
(B^), which depict the mixed bosonic fields of photons 
and collective atomic excitations. The eigen-frequencies 
of the polaritons A and B are 



a;2)2 + 16g2a;oW,(23) 



16.g2a;ow.(24) 



It is straightforward to see that uj\ < when g > gc ^ 
That is, the eigen-frequency uja of mode A be- 
comes a complex number, which means Hamiltonian (j22p 
will be non-Hermitian in the coupling region oi g > gc- 

The relations between the operators {a, 6, a\ b^} and 
{A, B, A\B^ are given by 



a' 
6t 



fiA^- 
hiA^ 



f2A- 
h2A 



hB^ - 
^h^B^ 



f4B, 
'hiB, 



(25) 



where the concrete forms of coefficients fi and hi (i — 
1, 2, 3, 4) have been given by Ref. Here we only give 
the detailed forms of ft in the Appendix. 

From Eq. , we can see that the ground state of the 
photon-dressed atomic ensemble in the polariton repre- 
sentation is |G) = |0)yi (8) |0)b = \00) AB- Making use 



of Eqs. ([20)1 and ([25)) . we can obtain the photon number 
variance 



7 = 2/2/1 + 2/1/1 + (/1/4 + f2hr 



(26) 



In the normal phase, all the coefficients fi {i = 1, 2, 3, 4) 
are real, then the photon variance is a positive number, 
which implies the coherence of the extra atom will vanish 
with time. 

We have mentioned that Hamiltonian (PT|) of two- 
coupled harmonic oscillators can not work well in the 
over-strong coupling region (g > gc). This is because 
the approximation -^/l — b'^b/N k. 1 for the original one 
[Eq. ([T^] can not make sense in this region. Thus, we 
need to consider a different approximation for Eq. (jl2p 
when g > gc- 



B. Dynamic sensitivity in super-radiant phase 

Physically, when the atom-light coupling becomes 
stronger and stronger, the coupled system will acquire 
a macroscopic excitations of atomic ensemble. And then 
the system enters into a super-radiant phase when g > gc- 
In this situation, the low-excitation approximation is no 
longer valid. We can use the coherent state |/3) of the col- 
lective atomic operator b to depict these kinds of macro- 
scopic excitations [1]. To achieve the effective Hamilto- 
nian over such background of macroscopic excitations, we 
need to do the displacement [1, Q 



(27) 



(or alternatively, 6^ — 6'^ -I- \^)- Correspondingly, we 
also displace the optical field by 



(28) 

\/a). Here a'^ and b'^ 



(or alternatively, 
describe quantum fluctuations about the semiclassical 
steady state [3; elsewhere, ^/a and y/P describe the 
macroscopic mean fields above gc in the order of 0{'\fN) 
5]. Then Hamiltonian (IT^ becomes 



Hg 



uJo W^b' - v/^(6't +b')+p 
+UJ [a''' a' + \/a(a'' + a) + o 



where 



{b'^Vi+V^b'-2^^), (29) 



VC=^Jl- [d^d - V^(dt + d)]/iN - P) 



is introduced. In the thermodynamic limit N ^ 00, for 
Eq. (|29p . we follow Emary and Brandes (Bj: expand the 
square root ^/J and keep terms up to the order of N'^ in 
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the Hamiltonian. Then through choosing the appropriate 
displacements 



with /i = LOLxJa/Ag^^ we can diagonahze Hamihonian ([2^ 
as 



by the Bogoliubov transformation 



= /(A't + f!,A' + f^B'^ + /IB', 



(30) 



(31) 



where the coefficients f- and /i- (i = 1, 2, 3, 4) have been 
given in Ref. Q • Here we only give the detailed forms of 
in the Appendix. 

The eigen-frequencies and lo'^ of the polaritons de- 
scribed by the operators A' and B' are given by 




It is known that if the coupling strength g exceeds the 
critical value gc, both the above eigen-frequencies are 
real, but not in the region of g < gc- Namely, when 
g > gc, Hamiltonian pop is Hermitian. 

In the super-radiant phase, the ground state \G) — 
\00)a'B' satisfies A'\G) = B'\G) = 0. Similar to the nor- 
mal phase, we can calculate the photon number variance 
in the super-radiate phase as 



7 = 2/f 



r, fl2 r/2 
J 4 



+ (/{/i + /2/3)' 



(34) 



Compared with the case of normal phase, the displace- 
ment a of the photon operator appears in the photon 
number variance. 



IV. PHOTON NUMBER VARIANCE FOR 
LOSCHMIDT ECHO 

We have separately calculated the LE of the photon- 
dressed atom ensemble perturbed by an extra atom in 
two quantum phases: normal phase and super-radiant 
phase. Our calculations are based on the short time ap- 
proximation, but it can cover the main character of the 
QPT of the photon-dressed atomic ensemble induced by 
the extra atom. As follows, we illustrate the LE versus 
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FIG. 2: (Color online) 3D diagram of the LE plotted as a 
function of the time t and the coupling strength g both in 
the normal phase (the left panel) and in the super-radiant 
phase (the right panel). Here, in unit of lj, ljo = L44ci;, 
S = Qs/As = O.OOlo; {As = O.lcj, gs = O.Olui), the critical 
point Qc = yJiOLoajl ~ 0.6aj, the number of atoms N — 100. 



the coupling strength g and time t by plotting its three- 
dimensional (3D) contour. 

Figure [5] shows the LE as a function of the time t and 
the coupling strength g in the normal and super-radiant 
phases. It is obvious that the LE, which is calculated 
from Eqs. (fT9|) . (f26|) . and (p4|) . will have a sudden change 
near the critical point. Its decay is highly enhanced at 
the critical value gc- In the normal phase, the LE decays 
rapidly to zero as the enlarged coupling strength g of the 
photon-dressed atomic ensemble approaches the critical 
point gc- In the super-radiant phase, similarly, the LE 
decays faster as the parameter g decreases to the critical 
point gc- Then the coherence of the extra atom is very 
sensitive to the dynamical perturbation of the photon- 
dressed atomic ensemble near the critical point. 

Meanwhile, in the vicinity of the critical point, the 
coherence of the extra atom decreases to zero sharply 
with time at fixed point of g- The more nearly the work 
point g approaches the critical point g^ the sharper the 
decay of the decoherence of the extra atom is. During 
this process, the detected atom evolves from a pure state 
to a mixed one. Therefore, we can measure the QPT 
of the photon-dressed atomic ensemble by exploring the 
coherence of the detected atom in the photon-dressed 
atomic ensemble. 

Figure [3] shows the LE at a fixed time {cot = 100) for 
the photon-dressed atomic ensemble in both the normal 
and super-radiant phases. Contrary to the case of the 
transverse field Ising model, the LE in the present sys- 
tem will not approach 1 when the coupling strength is 
much more than the critical point (seen from Fig. [3]). 
The reason is that a large displacement i/a a gVN ap- 
pears in the super-radiant phase and will increase as the 
coupling strength increases. That means a small dispar- 
ity (i5a^a) in the initial Hamiltonian in the super-radiant 
phase may lead to a large difference (e.g., the decoher- 
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FIG. 3: (Color online) The cross section of the 3D surface of 
the LE in Fig. [2]at ujt — 100. For other parameters see Fig. (2] 

ence factor will decay faster) after period of long-enough 
time. As pointed out in Ref. the so-called quantum 
chaos always appears in the super-radiant phase. 

It follows from Eqs. (dH), and jH]) that, the LE 
is independent of N in the normal phase. However, the 
LE depends on the number of the atoms N in the super- 
radiant phase via ^/a oc ^/N. In Fig.HJ the LE is plotted 
as a function of the coupling strength g with N — 100, 
1000, and 10000 respectively. It can be observed from 
Fig. |4] that the LE line decays faster and faster in the 
super-radiant phase as the atom number N increases. 
The reason is the same as that mentioned above. The 
photon number variance 7 proportional to the decay rate 
for the decoherence of the extra atom increases as N 
increases via approximately 

7 oc a oc g^N. (35) 

Accordingly, the LE decreases with the form 

InLoc-g^iV (36) 

in the super-radiant phase. Thus, as 00, the decay 

of the LE will be strongly enhanced at the critical point. 



V. ANALOG TO CLOUD CHAMBER 

Now we can address the similarity of sensitive dynam- 
ics between the present system and the classical cloud 
chamber. In classical cloud chamber, when a charged 
particle (or a dust) flies into the cloud chamber, which is 
filled with supersaturated and supercooled water or al- 
cohol, the water or alcohol vapor will condensate around 
the flying charged particle (or a dust) and form a liquid 
droplet, then a track is left. During this process, as a 
result of the sensitivity in response to the extra particle, 
the supersaturated vapor staying in the vicinity of the 
classical phase transition experiences a classical phase 
transition, transiting from vapor to liquid. 
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FIG. 4: (Color online) The LE of the systems for different N 
at Lut = 100. In normal phase, the LE is independent of N. In 
super-radiant phase, N — 100 (sohd line), 1000 (dashed hne), 
and 10000 (dotted line), respectively, from up to bottom. For 
other parameters see Fig. (2] 



In the present investigation, similarly, there exists very 
sensitive dynamics of the photon-dressed atomic ensem- 
ble when a far-off-resonant atom goes through the cav- 
ity. In view of the Stark effect, the far-off-resonant atom 
shifts the frequency of the cavity field. We assume that 
the photon-dressed atomic ensemble is initially prepared 
in a state near the quantum critical point of the QPT 
of the Dicke model. Then the frequency change induced 
by the far-off-resonant atom will lead the Dicke model to 
cross the quantum critical point, resulting in a sensitive 
dynamics of the LE. This quantum effect is similar to the 
classical phenomenon in the realistic cloud chamber that 
the vapor in the cloud chamber will condensate around 
the microscopic detected particle after experiencing the 
classical phase transition. Therefore, it is possible to re- 
alize the quantum version of the cloud chamber effect 
through observing the sensitive change in the LE of the 
photon-dressed atomic ensemble. 

Here, the enhancement of the decay of LE or its sudden 
change can be regarded as an indicator of the one-atom 
induced QPT to detect the passage of the atom. This 
fact properly resembles the cloud chamber effect. In this 
analogy, the photon-dressed atomic ensemble, which can 
be tuned to the vicinity of the QPT point, behaves as 
the supersaturated vapor in the classical cloud chamber, 
while the enhancement of the decay of LE just resembles 
the transition from vapor to liquid. 

Indeed, the LE in our paper is obtained from the de- 
coherence factor for time evolution of the extra atom, 
but it actually represents the "mark" of this atom on 
the "cloud chamber" — the photon-dressed atomic en- 
semble. An obvious reason is that the LE only depends 
on the parameters of the "chamber" and, thus, is an in- 
trinsic quantity of the chamber. Especially, the extra 
atom can only provide a small perturbation; thus, the 
LE is independent of the detected particle. In most of 
the references we cite, the LE can be defined without the 
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detected particle by the chamber. It is only in our own 
paper [2| where the detected particle is introduced and 
it is proved that the decoherence factor of the detected 
particle is just the LE of the chamber. Thus, the LE is 
obviously the mark of the detected particle left in the 
chamber. 



10775048, and the National Fundamental Research Pro- 
gram of China with Grants No. 2006CB921205 and No. 
2007CB925204. 



Appendix A: Coefficients of Bogoliubov 
transformation 



VI. CONCLUSION WITH A REMARK 



1. Normal phase 



In summary, based on the QPT of the Dicke model, 
we have proposed a quantum critical model to display 
the ultra-sensitivity of dynamic evolution of a QPT sys- 
tem of a photon-dressed atomic ensemble. We have also 
pointed out the analog of this one-atom induced QPT 
to the cloud chamber based on QPT. Frankly we have to 
point out that such a model can not be implemented eas- 
ily with the generic AMO system, since the two-photon 
term could not be simply ignored in the over-strong 
coupling limit 12]. However, our present study is still 
heuristic and the toy model covers the principle ideas 
for QPT inducing the cloud chamber-like effect. Fur- 
thermore, with the great development of solid quantum 
device physics, the Dicke model may be realized in some 
solid-state systems such as the super-conducting quan- 
tum circuits and the nano-mechanical resonators inte- 
grated with some qubit array systems. 

Finally, we would like to mention a reference [l^, in 
which an effective Dicke model was derived in a multilevel 
atomic ensemble. In this reference, the two-photon term 
may be safely ignored originally; thus, the modified 
Dicke model based on such a practical setup may be used 
to display the QPT phenomena we found in this paper. 



Acknowledgments 

We would like thank Shuo Yang for helpful discussions. 
This work was supported by the National Natural Science 
Foundation of China with Grants No. 10935010 and No. 



The coefficients of Bogoliubov transformation in the 
normal phase are 



1 cos 6* 

/l,2 = 77 , (UJ±UJAj, 



2 ^ujujA 
1 sin 



{lo±ujb), 



(Al) 



2 y/uJuB 

where the rotating angle in the coordinate-momentum 
representation 9 is given by 



tan20=^^^^. (A2) 



2. Super-radiant phase 

The coefScients of Bogoliubov transformation in the 
super-radiant phase are: 

f, 1 cos 6*' 

f 1 sin 61' 

where the analogous rotating angle 9' is 



tan 29' 



(A3) 



(A4) 



[1] S. Sachdev, Quantum Phase Transition (Cambridge Uni- 
versity Press, Cambridge, 1999). 
[2] H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. 

Sun, Phys. Rev. Lett. 96, 140604 (2006). 
[3] J. Zhang, X. Peng, N. Rajendran, and D. Suter, Phys. 

Rev. Lett. 100, 100501 (2008). 
[4] K. Hepp and E. H. Lieb, Ann. Phys. (N.Y.) 76, 360 

(1973); Phys. Rev. A 8, 2517 (1973); Y. K. Wang and 

F. T. Hioe, Phys. Rev. A 7, 831 (1973). 
[5] C. Emary and T. Brandes, Phys. Rev. Lett. 90, 044101 

(2003); Phys. Rev. E 67, 066203 (2003). 
[6] P. Zanardi and N. Paunkovic, Phys. Rev. E 74, 031123 

(2006). 

[7] D. Rossini, T. Calarco, V. Giovannetti, S. Montangero, 

and R. Fazio, Phys. Rev. A 75, 032333 (2007). 
[8] J. Zhang, F. M. Cucchietti, C. M. Chandrashekar, M. 



Laforest, C. A. Ryan, M. Ditty, A. Hubbard, J. K. Gam- 
ble, and R. Lafiamme, Phys. Rev. A 79, 012305 (2009). 
[9] L. C. Wang, X. L. Huang, and X. X. Yi, Phys. Lett. A 
368, 362 (2007). 

[10] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[11] Z. P. Karkuszewski, C. Jarzynski, and W. H. Zurek, 
Phys. Rev. Lett. 89, 170405 (2002); F. M. Cucchietti, 
D. A. R. Dalvit, J. P. Paz, and W. H. Zurek, ihid. 91, 
210403 (2003); R. A. Jalabert and H. M. Pastawski, ihid. 
86, 2490 (2001); T. Gorin, T. Prosen, T. H. Seligman, 
M. Znidaric, Phys. Rep. 435, 33 (2006). 

[12] K. Rzazewski, K. Wodkiewicz, and W. Zacowicz, Phys. 
Rev. Lett. 35, 432 (1975). 

[13] G. Liberti and R. L. Zaffino, Phys. Rev. A 70 , 033808 
(2004); Eur. Phys. J. B 44, 535 (2005); Y. Li, Z. D. 
Wang, and C. P. Sun, Phys. Rev. A 74, 023815 (2006); 



8 



D. Tolkunov and D. Solenov, Phys. Rev. B 75, 024402 
(2007); G. Chen, X. Wang, J. Q. Liang and Z. D. Wang, 
Phys. Rev. A 78, 023634 (2008); Y. Li and Z. D. Wang, 
[arXiv:0904.4730 
[14] H. Frohlich, Phys. Rev. 79, 845 (1950); Proc. R. Soc. 
London, Ser. A 215, 291 (1952); Adv. Phys. 3, 325 
(1954). 



[15] S. Nakajima, Adv. Phys. 4, 363 (1955). 

[16] T. Holstein, H. Primakoff, Phys. Rev. 58, 1098 (1940). 

[17] M. Wagner, Unitary Transformations in Solid State 

Physics (North- Holland, Amsterdam, 1986). 
[18] F. Dimer, B. Estienne, A. S. Parkins, and H. J. 

Carmichael, Phys. Rev. A 75, 013804 (2007). 



